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Abstract 

In this paper the several aspects of the symmetry in gauge theories at high tem- 
peratures are discussed. The metastable Zj^ bubbles in the SU{N) gauge theories with 
fermions may have, generically, unacceptable thermodynamic behavior. Their free energy 
F (X T"^ with a positive proportionality constant. This leads not only to negative pressure 
but also to negative specific heat and, more seriously, to negative entropy. We argue that 
although such domains are important in the Euclidean theory, they cannot be interpreted 
as physical domains in Minkowski space. The related problem is connected with the anal- 
ysis of the high-temperature limit of the confining phase. Using the two-dimensional QCD 
with adjoint fermions as a toy model we shall demonstrate that in the light fermion limit 
in this theory there is no breaking of the Zm symmetry in the high-temperature limit and 
thus there are no Z^ bubbles. 



* On leave of absence from ITEP, B.Cheremyshkinskaya 25, Moscow, 117259, Russia. 
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1. Introduction. 



The Z AT symmetry of pure Yang- Mills theories plays an important role in the study of 
their thermal properties It is known both from perturbative studies 0, ^ and from 
lattice simulations that at high temperature the Ztv symmetry is spontaneously broken 
and the the Euclidean theory has N degenerate vacua distinguished by different vacuum 
expectation values of Polyakov line |1[] 

< L >=< —tr P e' ^0 ^"'^^ > (1.1) 

In the presence of any matter which transforms as the fundamental representation of 
SU{N) (for example quarks in QCD or quarks, leptons and Higgs particles in the elec- 
troweak theory) this Zn symmetry is no longer present and all but one of these N de- 
generate vacua become either metastable or unstable (depending on N and the number 
of flavors). Following work on the computation of the interface tension between phases 
of different Z^r "vacua" in pure gluonic theories P it has recently been argued that 
these metastable vacua may lead to interesting cosmological consequences. 

However it has been found in that Z^v domains, generically, have unacceptable 
thermodynamic behavior, for example, their free energy F oc with a positive propor- 
tionality constant. This leads not only to negative pressure but also to negative specific 
heat and, more seriously, to negative entropy which means that something is definitely 
wrong in our understanding of the hot gauge theories structure. Using the abovemen- 
tioned thermodynamical arguments one must conclude that the Z^v -bubbles, i.e. different 
Ztv phases coexisiting in space, do not exist in models with metastable vacua. However 
the thermodynamical arguments do not forbid the existence of the different degenerate 
Z N phases in the theories with unbroken Z n symmetry. It is completely unclear what will 
be wrong with the theory if one adds some matter in fundamental representation ? How 
this matter can destroy the existence of (meta)stable states ? In some sense the situation 
would be much more clear if one can assume that these phases cannot coexist in the space 
even in the case of unbroken Zn symmetry. The problem of the Zjsf bubbles have been 
recently discussed in [||. 

Recently Smilga argued in a paper |^ that this point of view may be indeed correct 
and that " different Z jv thermal vacua of hot pure Yang-Mills theory distinguished in the 
standard approach by different values of Polyakov loop average correspond actually to one 
and the same physical state" . He presented different arguments supporting this statement 
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including the possible role of the infrared divergences in the calculation of the surface 
tension of the walls separating different Z^r phases as well as the difference between strong 
coupling lattice SU{N) gauge theories, where Z^v-bubbles exist indeed and the weak 
coupling continuum limit. 

Besides these general arguments he considered an interesting example of the 1 + 1 
dimensional hot QED - the Schwinger model - where the similar problem appears. Instead 
of Z]\f = 7!-i[SU{N)/Zn-] different vacuum states in a pure SU{N) gauge theory one has 
Z = 7ri[f/(l)] different states in the Schwinger model. However it was shown in |^ that 
in this case there are no domain wall solutions with finite surface tension. 

In this paper we shall consider another 1 + 1 dimensional model which, contrary to 
the Schwinger model, shares some common features with realistic 3 + 1 dimensional gauge 
theories. This is the 1 + 1 dimensional SU{N) gauge theory with Majorana fermions in 
the adjoint representation with the action 

Sadj = I d^xTr[ - -^F^^F^" + i^^D^,^ + m^^] (1.2) 

which obviously has Z^v symmetry. The light-cone quantisation of this theory was con- 



sidered in a large limit in |]10|. The spectrum consists of closed-string excitations. 



Contrary to the 't Hooft model [|rT| with fermions in the fundamental representation of 
SU{N) describing the open-string excitations with the only meson Regge trajectory, in 
this theory there is an infinite number of the closed-string Regge trajectories and the 



density of particle states increases exponentially with energy [O, [|T3| 



n(m) ~ m "exp(/5//m) (1-3) 

which means that there is the Hagedorn temperature Th = (3j]^ and the model undergoes 
a confinement - deconfinement transition at this temperature which must be the simplest 
analog of the real confinement-deconfinement transition in QCD. The numerical value of 



(3h ~ {0.7 — 0. 75) Jit/ (g'^N) in the large limit was calculated in [|T^]. The same picture 



was obtained in a recent paper [|1^ where the Hagedorn spectrum was obtained for a 1 + 1- 



dimensional QCD with adjoint scalar matter. The numerical value of inverse Hagedorn 
temperature in this case is Ph ~ (0.65 — 0.7) \Jtx /{g'^N). 

In a string theory language the Hagedorn transition occurs due to the fact that some 
winding modes in the imaginary time direction become tachyonic at high temperature 



[jT5[, |T6|) (see also IT^l)- In a string description of the gauge tyheory this winding modes 
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are generated by the Polyakov line operators 



1 fkP 

Lk{x) = — tr P exp (z / Ao{t, x)dr) (1.4) 
iV JO 



wrapping k times around the imaginary time direction r. In a very interesting paper |TS 
Polchinski studied the high-temperature hmit of the confinig phase and calculated in large 
N limit the mass of the tachyonic winding modes at high tempreatures. His method was 
used by Kutasov in |[T^ to study the stabilty of the confining phase in the two-dimensional 
QCD coupled to adjoint matter. We would like to note that both Polchinski and Kutasov 
analysis were based on the form of effective potential which leads to the existense of the 
Zj^ bubbles and these two problems - bubbles and tachyonic winding modes seem to 
be ultimately related. 

In section 2 we shall discuss the effective potential for Polyakov line and the thermo- 
dynamical properties of the Zn phases. In section 3 we consider the Polchinski method 
and discuss its relation to the existense of the Z^v phases. In section 4 the 1 + 1 gauge 
theory with adjoint fermions will be considered. Using the results obtained in [IS] we shall 
demonstrate that there are no coexisting Zn phases in this model by the same reasons 
which have been found by Smilga for Schwinger model j^. Moreover, there is no Z^r 
symmetry breaking in this model in te light fermion limit m << y/g'^N. In conclusion we 
shall discuss the obtained results and unsolved problems. In particular we shall consider 
the possibility that the metastable Z^ phases (if they do exist in the four-dimensional 
theories) may have interpretation as states with inverse population, i.e. with negative 
temperatures. 



2. Zjy domains in gauge theories. 



Let us briefly review the origin of the Z^ structure in gauge theories. It is most 
illuminating to begin with the Hamiltonian theory in = gauge and impose Gauss' 
Law DiEi — gip^tp = as a constraint on the Hilbert space of states. The projection 
operator P onto these states is just the projection operator onto gauge invariant states: 

P^r^au^ (2.1) 
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where the integral is over all gauge transformations g{r) with g{r) e SU{N) using the 
Haar measure Vg. Ug is the representation of the gauge transformation g on the Hilbert 
space of states. The partition function at nonzero temperature is given by 

Z = Tr [e-^"p} = jvgTr {e'^^Ug} (2.2) 

If we compute this trace in a basis {l^j,^)} where ^ represents an appropriate fermionic 
state then 

Z^J VA,Vi j Vg (A, e|e-^^t/s| A, (2-3) 

One can now proceed with the usual derivation of the Euclidean functional integral at 
nonzero temperature except to note that the presence of the factor Ug modifies the bound- 
ary conditions on both the gauge fields and the fermions. Thus apart from an overall 
normalization 

j Vg{r) [ VAi{f,T)ViP{f,T)Vip\f,T)e- ^0 '^^ ! 
J Jb.c. 

B.C. : A(r, (3)^g [A{r, 0)] P) = -g 0)] (2.4) 

where L^; is the usual Euclidean Lagrangian for a gauge theory coupled to fermions, 
g[A] = gAg~^ — dgg^^ and glip] = gip. In other words we derive the usual functional 
integral in = gauge but where the periodic (antiperiodic) boundary conditions are 
modified to be periodic (antiperiodic) up to an arbitrary gauge transformation which we 
then integrate over. 

It is of course possible to remove these strange boundary conditions by performing 
a gauge transformation in the functional integral. For each value of the integrand g we 
may introduce any gauge transformation V{f,T) with the property that V{f,0) = I and 
V{f, P) = g^^{r). This will force the introduction of a temporal gauge field Aq = doV V~^. 
The integral over g will then become an integral over Aq with the appropriate measure. 
In fact if we integrate over all possible such V^s we recover the usual path integral over 
all gauge fields A^{f,T) with periodic boundary conditions and over all fermionic fields 
with antiperiodic boundary conditions. 

If we consider spatially constant gauge transformation g G Z^v i.e. g = exp{27iik/N) x 
/ (where k is an integer and I is the identity matrix). Then g[Ai] = Ai but g[^] = 
exp{2TTik/N)ijj. Thus in the absence of fermions there are N degenerate vacua. When 
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fermions are present, however, a value of g E corresponds to a path integral in which 
the fermions have the "twisted" boundary conditions = — exp{2TTik/N)ilj{0). One 
can say that this describes fermions having an imaginary chemical potential (for more 
detailed discussion see, for example 



The main tool for analyzing the Zn structure of gauge theories is the calculation of 
the effective potential in a constant, background temporal gauge field Aq. For definiteness 
we begin by considering an four dimensional SU{N) gauge theory with Nf flavors . of 
massless Dirac fermions. Although at finite temperature it is not possible to choose a 
gauge in which = it is possible to choose a gauge in which Aq is independent of 
(Euclidean) time r and in which it is diagonal. Aq can then be written as 



<d = /3AQ 



\ 



\ 



On J 



(2.5) 



where for SU (N) one has 9i + . . . + 9n = (mod In / f3) and there are — 1 independent 
di - the number of independent diagonal generators of SU{N) (i.e. the elements of the 
Cartan subalgebra of the Lie Algebra of SU{N)). In this gauge the Polyakov line is given 
by 



1 ^ 
-Ye'' 



[2.6) 



The effective potential for Aq has been calculated up to two loops [|], ^ Here we 
shall consider only the one loop result P, H]. For gluons the effective potential is 



Vn 



and for each fermion flavor in fundamental representation it is 




i-, - ■ - On 








/ mod 2 


^/mod2. 



-N - 



45 



12 



1=1 























/ mod 2 J 





(2.7) 



(2.8) 



The values of the effective potentials at zero field Vg(0, . . . , 0) = — (7r^T'^/45)(A^^ — 1) 
and Vf(0, . . . , 0) = — (7/4)(7r^T'^/45)A^. It is easy to see that it is the free energy density 
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of an ideal gas of gluons and fermions at a temperature T which is equal to — tt^T^k/QO, 
where each bosonic degree of freedom contributes 1 to k and each fermionic degree of 
freedom contributes 7/8 to k. Thus the gauge fields contribute 2{N'^ — 1) and the fermions 
contribute 4:{7/8)NfN to k which reproduces ^^(0, • • • , 0) and Vf(0, . . . , 0). 

The gluon effective potential (|2.7|) has symmetry 



2nk 



(2.9) 



where integer k = 0, . . . , N — 1 must be the same for all 6i. Then = (mod27r//?) 
and without fermions there are N minima at 6i = {27r/f3N)k where k = 0, . . . , N — 1. 
The fermion potential ( ^.81 ) obviously violates the Z^v symmetry ( [^.9D and only 6'j = is 
the global minimum of the total effective potential. All other Z^v vacuua become either 
local minima (metastable states) or unstable states depending on the number of fermion 
flavors Nf. 

Furthermore it will be sufficient for us to consider the case in which 6 has the form 



n 



/I 



V 



:2.io) 



-{N-1)J 



The total free energy as a function of q is the sum of a gluon ( p.7| ) and Nf fermion 
effective potentials and can be written as 



F{q) 



-iN-l){VG{q) + NfVF{q)) 



1 

45 



1) + ^-NNj 



TT^T^ 3 ' 

where the bosonic and fermionic contributions are given in terms of the function 



(2.11) 



I, "''mod 1> 



•^mod l) 



:2.i2) 



by 



Vciq) = /(g), Vp{q) 



16A^-1 •'^A^ 2' 



A^- r ^A^ 2 



q) (2.13) 



Notice that the function / is periodic with period 1. Thus Vciq) is periodic with period 
1 but Vpiq) is periodic with period A^. Note also that a constant has been added to Vp 
so that it vanishes at g = which is the perturbative vacuum of the theory. The last 
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term in (|2.11 ) is the free energy density of an ideal gas of gluons and Nf fermions at a 
temperature T. 

Now let us repeat the arguments which were used in to demonstrate very serious 
problems arising in the thermodynamical description of the metastable Z^- vacuua. For 
example for = 3 it is easy to show that the metastable minimum at g = 1 remains 
metastable for Nf < 18 at which point it becomes unstable. Notice, however that for 
Nf > 3 the free energy density becomes positive in these metastable states. Since 
F ocT^ this poses a very serious problem which we shall now discuss. 

First note that the positivity of the free energy at nonzero values of q is entirely due 
to the fermions. For integer q < N/2 one gets 



F{q integer < N/2) 



NNf^td---:^ 



180 



AT^-l 
45 



[2.U) 



This will be positive provided 

Nf> 



45 N 




1 - 



2f 

N^ 



7 

180 



-1 



(2.15) 



For = 3 and q = 1, for example, we find that F is positive (and propotinal to T^) for 
Nf > 3. For even iV and g = N/2 one has positive F{N/2) for Nf > {N^ - l)/2iV. This 
point will be a local minimum for Nf < N. Thus for = 4 and N2 = 2, 3 the point 
g = 2 is both a minimum of F and has a positive value F propotional to T^. In more 
details this was discussed in ||^. 

It follows from the above discussion that for a large class of models with both gauge 
fields and fermions we can write the free energy density as 



[2.16) 



Such a situation is impossible for the metastable state of a real physical system. To 
see this let us remember that the free energy of any physical system at temperature T is 
defined as 

F(T) = -Tln^e-^"/^ (2.17) 

n 

Shifting En by a constant C one can add C to F{T), but the T dependent part must be 
negative as one can see defining energy levels En in a such way that the ground state energy 
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Eq = and F{T) = — Tln(l + ...)< 0. In our case we got positive F{T) which leads to 
physically senseless thermodynamic quantities, namely the negative entropy density 

^ E-F dF(T) , , , , , 

S = = = -4|7|T^ (2.18) 

the negative internal energy density 



E = F + = -3|7|T^ (2.19) 



the negative pressure 



P = -W (2.20) 

and the negative specific heat 

c=-12|7|r^ (2.21) 

Such a metastable vacuum thus has not only a negative pressure but also negative specific 
heat and worst of all a negative entropy. It is clear that no physical systems with positive 
temperature of this type can exist. However one cannot exclude the metastable states 
with inverse population (the well known examples are lasers). We shall briefiy discuss 
such an intriguing possibility in a conclusion. 

There are interesting cases in which the free energy density at the metastable minimum 
has the correct, negative sign for F . One such example is given in Ref. in which the 
standard electroweak model is considered well above the QCD phase transition point. 
In this case the base free energy density of the leptons, the Higgs and the weak gauge 
bosons contribute to the total free energy density and make it negative. It is clear, 
however, that if a subsystem of the full system (namely the quarks and gluons) has the 
disease discussed above, namely a positive free energy density which grows like T^, then 
including the leptons, Higgs and weak gauge bosons which couple weakly to it cannot save 
the situation. This will be clarified below but if we imagine a system whose total entropy 
is positive but that some identifiable subsystem has negative entopy then any statistical 
description of the full system fails since the subsytem has no states available to it. More 
discussion about thermodynamical properties of these states as well as difficulties arising 
in interpretation of domains in Minkowski space can be found in 0. 

Let us note that the real problems arise only when we are assuming that this metastable 
states appear in our space as Zat bubbles. To get such bubbles it was usually assumed that 
one can use the effective potentials ( |2.7| ) and ( |2.8| ) not only for constant Aq for which the 
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potentials had been really found in |Q] and but also for spatially dependent Aq{x). Let 
us note that symmetry (neglecting the matter in a fundamental representation) exists 
only for coordinate independent part Aq and generally speaking it is not necessary that 
potential for nonconstant modes Ao{x) is the same. Moreover, as it has been demonstarted 
by Smilga in Schwinger model |^ and as we shall demonstrate later in a 1 + 1 QCD with 
adjoint matter, there are situations when the total effective potential is the sum of two 
independent potentials for constant and nonconstant modes. In this case the zero mode 
Ao becomes quantum mechanical variable and no Zn bubbles exist in space - however 
the price for this is the unbroken Z^ symmetry. Before we shall start this discussion let 
us consider how one can use the effective potential ( pT7| ) to analyse the high-temperature 
limit of the confining phase following the ideas suggested by Polchinski in [113] . 



3. Z]y phases and high-temperature Umit of the confinig phase. 

Let us consider the two-point correlation function of Polyakov lines (|1.4|) at low tem- 
perature in a pure gluodynamics 

(Lfc(f)L_fc(0)) ~exp(-Mfc(/5)x), a; ^ oo (3.1) 

The correlation function vanishes at infinity because of the confinement, so the expecta- 
tion value of Polyakov line is zero (L^,) = and Zj^ symmetry is unbroken. The usual 
interpretation of the Polyakov line is the world line of an external source in fundamental 
(for /c = 1) or in higher [k > 1) representations. However there is a dual description when 
one can consider compact Eucledian time r as a spatial coordinate and Lk is the creation 
operator of a winding state with a electric flux in the periodic direction and Mk{l3) is a 
temperature dependent mass of the winding state. In the deconfinement phase the wind- 
ing modes become tachyonic and the theory makes a transition to a phase with broken 
Zn symmetry. 

Polchinski in |]18| used the high-temperature effetctive potential (p77|) to calculate the 
mass of these tachyonic states. To find them he considered the effective action 

1^2 2 
24^2^3 ~ ^fc)mod 2^ [27^ " (^i " ^fc)mod 2^] (3-2) 
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The gradient term corresponds to the square of the electric field E'^ in the bare action 
and the second term can be obtained from the effective potential ( p.7|) after omiting 9 
independent first term. In the large N limit one can introduce the normalized density 



1 



(3.3) 



In the N oo limit symmetry is transformed into U{1) symmetry 9^9 + const. It 
is easy to see that L^^x) are the Fourier coefficients of the density p{9, x) (see ( p.6|) ) 



Zn 



Lkix) 



1 + ^ Lk{x) exp{—ik9) + L_k{x) exp{ik9) 

k=l 

ike - 1 — , _ , , , i .11- 







1 " 1 rkp 

p{9, x)e"'^d9 = — y exp {ik9i{x)) = — trPexp {i / Ao{t, x)dT) (3.4) 
N ~[ N Jo 



The action ( p.2|) takes the form 



eff 



N 
2^ 



(fx 



2lT 



d9 



p{9, X 



id^'Vp{9,x)y 



247r2/53 



/ d^x r d9^ r d92p{9i,x)p{92,x){9i - 92)[2tt - {9^ - ^2)]' 
J Jo Jo 



(3.5) 



where dg^ is defined as dg^ exp(ik9) = exp(ik9)/ik and in large limit we must keep 
g'^N fixed. The minimum of the potential is when all the eigenvalues are equal which 
gives us the spectral density in high-temperature phase with broken (here is U{1) in 
large N limit) symmetry 

Pbrokcn(6', x) =5(9- 9o) (3.6) 

for some 9o- The symmetric confining phase is defined by the U{1) invariant distribution 

(3.7) 



pc[t^,X) 



1 

2^ 



which is unstable. One can easily find that in quadratic approximation in the action 
(13.51) takes the form 



oo „ 



X 



k=—oo ' 



2g'^Nl3k 



247r2/33 



(3.^ 



where Vk is the Fourier transform of the potential V{9i — ^2) 

l2 



9. 



2 /mod 2it\ 



1 



27r Jo 



Vk = — I d9e'''^V{9) = — I d9e'^''9'\2n 



2tx Jo 



2-K 



AkQrai 



h-02)L, 2.[27r-(^^i- 

24 , , 

= (3.9) 
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and one gets the tachyonic winding modes with masses 

Ml = -2g^N/7r^(3^e (3.10) 

One can make an assumption that the same spectrum of the tachyon masses must be 
reproduced in the string theory (if any) describing QCD. Let us note that this spectrum 
is different from the spectrum obtained in the usual string theory o^'M^ = —n/6 + 
/3^/47r^a' where n is some effective constant proportional to the number of the world- 
sheet degrees of freedom (24 for critical bosonic string). Using ( |3.10| ) Polchinski made 
a conclusion that for QCD string the effective number of degrees of freedom grows with 
temperature as neff{j3) ~ g'^{l3)N / {3"^ and the main conclusion of his analysis is that a 
string theory describing QCD in the large N limit must have a number of world-sheet 
degres of freedom which diverges at short distances [|l^] . 

Let us note that this analysis can be repeated even in the case when there are Nf 
fermions in the fundamental representation. Then using the ( f^.8[ ) one can see that the 
fermion contribution to the effective action (|3.2| ) will be 

Of - 



127r2/33 



J d'x E - ( i^i + ^)mod2. - vr)'] ' (3.11 

i=l 



This term corresponds to the linear in p{6, x) contribution to the effective action for the 
density p 

Sf = / d'x I^d9pi9,x)i9^ - = ^^f^j '^"'^^M^) (3-12) 

Including this term into ( p.8|) we get 



CO „ 

k=—oo 



1 - - 1 Nf (-)'= 



(3.13) 



and it is easy to see that linear terms do not change the values of the tachyon masses 
(|3.10|) , but shift the fields 

L,^L,+ ^-J-^ (3.14) 

For small Nf/N this shift is small and the quadratic approximation for the total action 
( |3.5| ) is still valid. One cannot use this approximation when Nf ^ N, i.e. precisely when 
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the suppression factor for nonplanar diagrams Nf/N is of order one and nonplanar dia- 
grams have the same order of magnitude as the planar ones. In latter case the connection 
between large N gauge theory and a string theory is much more problematic. 

Thus we see that effective high-temperature theory gives us an important informa- 
tion about (possible) string description of the low-temperature confinement phase. This 
information was based on the form of the effective potential which may lead in some 
situation to physically senseless results. So we must be very accurate in dealing with this 
potential and for this reason let us study the simplest model where one can hope to have 



confinement-deconfinement transition -1 + 1 dimensional QCD with adjoint matter |T0 

- nn. 



4. Two-dimensional QCD coupled to adjoint matter at high 
temperatures 



Let us repeat following Kutasov [|l^] the analysis of the previous section in the case of 
the 1 + 1-dimensional QCD interacting with Majorana fermions in the adjoint represen- 
tation described by the action 



S, 



adj 



dx I drTr 





(4.1) 



defined on a Eucledian space-time with periodic Eucledian time r ~ r + /? with periodic 
boundary conditions for gauge fields and antiperiodic for fermions. One can again choose 
the gauge when Aq is diagonal and independent on time r: Aq = (1//3) diag(6'i, . . . 9^)- 
Then the one-loop effective action for the 6i takes the form 



S, 



1 f ^ 



da 



i\2 



dx 



+ V{9^,...9 



N 



(4.2) 



where the effective potential is nothing but the determinant of the Dirac operator in 
adjoint reperesentation j^Dlf' in an external field 



V{9^ 



9 



N 



— In det 



N 



m 



J2 lndet[^^D^[9,-9j]+m] 



(4.3) 
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Using the proper time representation of the fermion determinant one gets 



^(^« - = dxYlH'l 7?exp(-^ - mV) cosfc(^. - 0,) 

yi0^ - ^i) = 77-3 f dx[{e + 7r)^od2. - vr]2, m = (4.4) 

Znp J 

and in the high-temperature hmit j3 —>■ one can neglect the mass term in the leading 
approximation. The sum over k is the sum over winding of the particle trajectory 
around the compact imaginary time and in the high-temperature limit /3 — one can 
neglect the mass term [| in the leading 1/(3 approximation. This representation of the 
determinant can be easily generalised to higher dimensions where one must substitute 
(ir/r^ by dr/r^^^ in the case of d + 1 dimensional theory. Also in the case of bosonic 
degrees of freedom the factor {—)^ is missing. One can easily reproduce (|2.7|) and ( p.8|) 
using the proper time representation of the boson and fermion determinants. 

Thus for small (3 the effective potential for constant 6i is known and assuming, as 
usual, that for slowly varying 9 one can simply substitute 9i by 9i{x) one gets the effective 
action (see [P^) 

1 „ Af jn Or N oo -I 

s^jf = ^ / + ^ / E(-) V - e^)] (4.5) 

Obviously this action has symmetry 9i ^ 9i + 27rm/A^, m = 1,2, ... N and we get 
the same Zjy bubbles as in four-dimensional case. One can also consider the density (|3.3|) 



p{9, x) and using the same method as in |T^ Kutasov obtained the effective action for 
Lfc(x) in quadratic approximation [| 



oo „ 

S = N^ J dx 

k=—oo 



(4.6) 



and got the masses of the winding modes 



M,^(/3-0) = ^(-)^ (4.7) 

7T 

where only the odd k winding modes are tachyonic. 

*it is easy to see that the corrections wih be by order of mf3 

^ To get the effective action in quadratic approximation one must simply repeat the same 
procedure as in a 3 + 1 dimensional case. Let us note also that because we need Fourier 
transformed effective potential Vk it is not necessary to take the sum over k in ([4.5|) to get V{9). 
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However in the massless limit the form of the effective potential (4^) which leads both 
to the existence of the domain walls (and in the case of the fermions in the fundamental 
representation they will be metastable bubbles) and gives the imaginary mass of the 
unstable winding states is wrong. As we shall demonstrate now in the massless theory 
there are no domain walls and one cannot make simple predictions about the condensation 
of the winding states - and all this because the potential for the x dependent modes of the 
fields di is not periodic contrary to the constant mode part which was really calculated in 
The two-dimensional fermion determinant on cylinder in arbitrary gauge field was 



calculated in [|T9| where it was shown that it is factorized into product of two independent 
terms. The first one is the periodical potential but for constant modes 9i, not for the fields 
6i{x). The second term depends on x-dependent fields 6i{x), however there is no reason 
why this part must have Z^v symmetry - and it has not. In result any configuration with 
domain wall has an energy proportional to the one- dimensional system length / dx and 
these configurations are absolutely irrelevant. 

To make this statement more clear let us repeat the calculation of the determinat of 
the two-dimensional Dirac operator and demonstrate the factorization property following 
the paper ||19|. We shall not consider here the case of an arbitrary gauge field A^, where 
Afj_ is a Hermitian matrix, instead we shall consider interesting for us case when this 
matrix is diagonal (see (|2.5|)) 



A^(r, x) 



Then the fermion action from (|4.1j ) takes the form 



fAl{T,x) 



\ 



A^^{T,x)j 



Sf = J dx drTr {i-ff,d^ + m) * + ^7^[A^^ 

J dx f dr J2 {^1,9^ + - ^j] + ^) ^ 



(4.9) 



and we see that the total determinant for the adjoint fermions 

deti^,D^,,^iA) +m) = Y[deti^,[td'^ + Af^] + m) 



(4.10) 



is the product of abelian determinants of fermions interacting with abelian fields A'^^ 
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As we discussed before in a high-temperature hmit one can omit the mass term and 
what we are trying to calculate now is the determinant of the Dirac operator for the 
massless fermion in an abelian gaauge field on the cylinder x or on the torus 

X if the space is a circle too. Let us consider the Hodge decomposition of the vector 
potential (in the sector with zero topological charge) 

At" = ^e^' + e^'^d^X + d"^ (4.11) 

where the last term is the pure gauge and can be neglected. The i?^ factors are the radii 
of two 5*1 and Rr = 13. In the case of the noncompact space, i.e. when we are considering 
the cylinder x R^ one can put R^ oo. Constant modes 9^ are the angle variables 
with periodicity 2ti and x is a coordinate dependent mode. Substituting the Hodge 
decomposition into Dirac operator and using the well known property of two-dimensional 
7 matrices 7^e'^'^ = i7'^75 one can get 

^7^(9^^ - iA^") = Z7^ (^d^' - -^^^ - k'^^a^x^ = ^e^^^7^ - -^^^^ e^^^ (4.12) 

Now let us consider a family of operators 

Vr = ie^^'^'-ff. Id" - -^6 A e"^^^ (4.13) 



We are looking for Di and Dq = 7^ [d'^'' — -^O^j is the operator in a constant field which 
determinant is periodic in 6'^. 



A very elegant formula was obtained by Blau, Visser and Wipf in |19] for a family 



of general first order elliptic self-adjoint operators depending on a parameter r of the 



type Or = y (7o/5't exp(r/''')Co exp(r/) were f{x) can be in general some matrix-valued 
function and Qr is the parameter dependent metric on the manifold. In our case (|4.13| ) we 
have f{x) = j^xi^) and fiat metric independent on r. We shall present here the result for 
this case only , more general expressions including general nonabelian field can be found 
in original paper [|19|. For a family of operators (|4.13| ) one gets 



-^IndeWr = - [ (fxxd\ (4.14) 
dr TT J 



^ Here we neglect the possible contributions of the zero modes. If they are present the 
determinant itself is zero and one has the expression for dei'V^/detKr, where is the matix of 
scalar products of zero modes. However if we consider the sector with zero topologicaal charge 
there are no zero modes and one can simplify life a little 
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which can be easily integrated and finally we get 

det^^id^ - zA^) = det^^ (^d^ - -^r^ exp (^i- J d''xxd\^ (4.15) 

and we see that the contributions of constant modes 6^ and nonconstant modes x{^) ^.re 
factorizable ! Thus the effective potential is the sum of two terms - one is the effective 
potential for the constant field which we had calculated and which is periodic in 6 with 
period 2tt. Another one is the potential for x-dependent part of the gauge field which is 
not periodic at all - this is the ordinary Schwinger mass term proportional to / d'^xxd'^x = 
il/A)fd^xe^,F^,il/d^)e^,F^,. 

Now let us consider the effective action when only Aq component in ( [4.8|) is nonzero 
and Aq depends only on x. This is precisely the case which is relevant both for studying Z^r 
bubbles and instability of the confining phase at high temperatures as we have discussed 
before. Then one can immediately write the effective action in which the constant modes 
9i are completely independent from x dependent fields Ai, where "means that we extracted 
the zero mode from A- fx). 



S, 



eff 



^ - + ^)mod 2^ - T^f + 

N 



i=l 



27r 



(4.16) 



where L is the size of the one-dimensional system. The symmetry of this action 
Oi ^ Oi + 2i\mlN , m = 1,2, . . . N does not affect the x dependent fields Ai{x) at all. 

It is clear now that for such potential any Zn bubbles will have infinite energy in the 
thermodynamical limit. To see it let us consider the field (see (|2.5| )) 

/I 

271 



Ao{x) 



(3N 



q{x) 



\ 



V 



(4.17) 



-(iV-1)/ 

where g = and q = 1 corresponds to two different vacuua. The effective action for 
one interpolating field q{x) is the following 



S[q] = {N 



2Tf_ 



27r 
J 



27rL 



mod 1 



■(4.18) 
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where q{x) = qo + qix) and zero mode go = (1/-^) Jq dxq{x). Let us consider the 
bubble, i.e. the region with q{x) = 1 with size R. The rest part of the space, i.e. region 
with size L — R has q{x) = and we assumed that both L and R are large in comparison 
with domain wall width 1/y/Ng'^, so one can neglect the contributions from the regions 
where q{x) interpolates between and 1. It is easy to see that in this case go = R/L and 
g(x) = q{x) — R/L and one must have R/L < 1/2, in opposite case it is better to say 
that there is a bubble of g = phase in the space with g = 1. Then the second and third 
term in S[q] contribute to the action 



27r 

S{L,R) = {N-l)j 



R{i - + {L - R){^r + Li^r 



{N-l)jR (4.19) 



and we see that the action is proportional to the bubble size R. To consider the single 
domain wall one must put R L sending the "anti"-wall far apart, then in the thermo- 
dynamical limit L oo and single domain wall is infinitely heavy. 

Thus there are no bubbles in this theory, moreover we are not allowed to have 
any Zjv vacua. The reason is very simple - in this theory the Zjv symmetry can not be 
broken spontaneously. The reason is obvious - due to the decoupling of constant modes 
from all other modes we have to average over all Z^ vacua and cannot restrict ourselves 
for only one - as it would be possible in the case without factorization. Let us consider 
for example SU{2) theory with the symmetry of the center Z2. The Polyakov line in this 
case is |^ 

and it is evident that < L{x) >= because averaging over go is factorized and 

< e^^""" >g,= (4.21) 

However the zero-mode go factors are cancelled in the two-point correlation function 

< L{x)L{0) >=^< cos(g(x) - g(0)) >= le<«»«'(°)> (4.22) 

which is non-zero and has finite limit at x 00. Thus we have < L{x)L{0) >^ 1/2 
when X —>■ 00 and at the same time < L{x) >= 0. The breakdown of the clusterization 
property is connected with the factorization of the zero mode go. This mode is absolutely 
delocalised and one can not use naive rule < A{x)B{0) >-^< A{x) >< B{0) > at a; — >• cxo. 

^Here we are considering the simplest case with the winding number k = 1. The generalization 
for arbitrary k is trivial 
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Because the symmetry is unbroken at high temperatures one can conclude that 
the analysis of the instability of the confining phase at high temperatures which was done 
in is incomplete. However we must remember that the factorization of the fermion 
determinant ( [4. 151 ) was obtained only for massless case. For massive fermions the effective 
action is 

oo 

Tr \n[^^D^ - m] = Tr In 7^^^ - Tr ^ (4.23) 

n=l WM-'^ ) 

and only in first term there is factorization, other terms may mix constant and variable 
modes. However their contribution is suppressed by powers (3m and one can think about 
their importance at temperatures T < m. It is still unclear what is the order parameter 
corresponding to the Hagedorn transition in the theory with light fermions with mass 
m < \J g'^N/n. In any case the usual picture of the Zn symmetry breaking is wrong in 
this case and one has the strange picture of unbroken Z^q with < L{x) >= 0, but with 
nonzero < L(oo)L(0) > which means deconfining. 
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5. Conclusion 



In this paper we tried to address some at first sight different but, from our point of 
view, related problems connected with the symmetry in the hot gauge theories. We 
demonstarted that allowing the nontrivial Zn vacua to exist one may have a real problem 
when studying the metastable states which have impossible thermodynamic properties. 
We tried to present arguments that such bubbles cannot exist as metastable states at all. 
However we did not explore in this paper another possibility, which seems unrealistic, 
but can not be excluded completely. This is the idea that if after all our attempts the 
metastable states will survive one must consider them as some kind of states with inverse 
population, i.e. the states with negative temperature! Let us remember that both 
negative specific heat and negative entropy proportional to and for negative T they 
will change the sign, i.e. in the case when they were negative at positive T they will be 
positive for negative T. This idea seems rather strange because we still have a paradox - 
what we have started from was the gas of quarks and gluons at high positive temperatures. 
How does one get the negative temperature is a big question. One may speculate that 
we can reach this region passing through infinte T. If such a thing would be possible 
one could imagine something like quark-gluon laser in hot gauge theories. Despite all it 
strangeness this idea deserves further analysis. 

We also considered connection between existence of Zn bubbles and the the insta- 
bilities of the confining phase at high temperatures. Using two-dimensional QCD with 
adjoint fermions as an example we had demonstrated that situation with the breaking 
of the Z]\f symmetry is not so simple as one could imagine when dealing with "simple" 
two-dimensional models. In fact, for massless case we have proved that the Z^ symme- 
try can not be broken at all, however the correaltion function of the two Polyakov lines 
< L{x)L{0) > is not going to zero when x oo and we are in the deconfining phase at 
high temperatures. It is unclear how the mass of the fermion affect this situation. We 
conjectured that the same physics must be correct for light enough fermions with mass 
m less then y/g'^N. For heavy fermions with > g'^N one may have another phase 
and may be even the breaking of global Z^. It is amusing that in this theory one has 
hidden supersymmetry precisely at = g'^N/n as was shown by Kutasov in [Q. Does 
it mean that SUSY point is a phase transition point and one has two phases in a QCD 
with adjoint fermions is a very intriguing possibility. 

It is also unclear how reliable is the effective action in the four-dimensional gauge 
theories at high temperatures. What is known is the effective action for the constant 
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field and one assumes that it is possible to substitute the constant field Aq by general 
x-dependent field Aq{x). We just saw, however, that this recipe is not good all the time - 
two-dimensional example is a good demonstartion of this fact. May be in four dimensions 
we also missed something ? In any case it seems that the problem of Zn phases and 
related with it the problem of the instabilities of the confining phase Q definitely deserve 
further investigations. 
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